Isospin asymmetric nuclear matter and properties of axisymmetric neutron stars 
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Pure hadronic compact stars, above a limiting value (~1.6 Mq) of their gravitational masses, to 
which predictions of most of other equations of state (EoSs) are restricted, can be reached from the 
equation of state (EoS) obtained using DDM3Y effective interaction. This effective interaction is 
found to be quite successful in providing unified description of elastic and inelastic scattering, various 
radioactivities and nuclear matter properties. We present a systematic study of the properties of pure 
hadronic compact stars. The /3-equilibrated neutron star matter using this EoS with a thin crust is 
able to describe highly-massive compact stars, such as PSR B1516+02B with a mass M=1.94j;Q jg 
Mo and PSR J0751+1807 with a mass M=2.1±0.2 M© to a la confidence level. 
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PACS numbers: 21.65.Cd, 21.65.Ef, 26.60.Kp 



The theoretical study of the nuclear equation of state 
is a field of research which ties together different areas of 
physics. Nuclear EoS is of great interest as its features 
control the stability of neutron star (NS), the evolution 
of the universe, supernova explosion, nucleosynthesis as 
well as central collisions of heavy nuclei. Extensive stud- 
ies in the past two decades of nuclear matter created 
at subnormal or supernormal density in heavy ion col- 
lisions have resulted in experimental constraints on the 
nuclear EoS of symmetric matter. Recent astrophysi- 
cal observations of massive neutron stars and heavy-ion 
data are confronted with our present understanding of 
the EoS of dense hadronic matter. We argue that the 
data from massive neutron stars and pulsars provide an 
important cross-check between high-density astrophysics 
and heavy-ion physics. The density dependence of nu- 
clear symmetry energy (NSE) obtained by using nuclear 
EoS plays an important role for modelling the structure 
of the neutron stars (NSs) and the dynamics of super- 
nova explosions since a series of observables (e.g. slope 
L of NSE, the value of NSE at nuclear density etc.) can 
be determined from the knowledge of symmetry energy. 
The stiffness of the high-density matter controls the max- 
imum mass of compact stars. New measurements of the 
properties of pulsars point towards large masses and cor- 
respondingly to a rather stiff EoS |1| characterized by 
symmetric nuclear matter (SNM) incompressibility 250- 
270 MeV or more. In a recent analysis of x-ray burster 
EXO 0748-676 (M= 2.10 ± 0.28 Mq) it is even claimed 
that soft nuclear EoS are ruled out . In addition, it is 
argued in Ref. that condensates and deconfined quarks 
Q may not exist in the cores of NSs. Recently, new ob- 
servations for the mass and the radius of compact stars 
have been obtained which provide stringent constraints 
on the EoS of strongly interacting matter at high densi- 
ties 0,11. 



We investigate the impacts of the compression modulus 
and symmetry energy of nuclear matter on the maximum 
mass of NSs in view of the recent constraints from the 
isospin diffusion in heavy-ion collisions at intermediate 
energies d, 01 • In the present work, the density depen- 
dent M3Y effective interaction (DDM3Y) Q which pro- 
vides a unified description of the elastic and the inelastic 
scattering [§, [H , various radioactivites [lll - fTsj and nu- 
clear matter properties [lB43: is employed to obtain 
EoS of the (3 equilibrated NS matter. This EoS is used 
to carry out a systematic study of the static as well as 
rotating NSs in view of the recent observations of the 
massive compact stars. 

As mentioned above, the nuclear matter EoS is calcu- 
lated using the isoscalar and the isovector components of 
M3Y interaction along with density dependence. The en- 
ergy per nucleon for isospin asymmetric nuclear matter 
is given by 



2l2 



3?i 

^ 10m 



F(X) + (^)(l-/3p") (1) 



* |royc.partha®gmail .coni] 
t abphy@caluniv.ac.ml 
■f .dnbQveccal.ernet.inl 



where is the isospin asymmetry with /5„, pp 

and p—pn+Pp being the neutron, proton and nucleon 
number densities respectively, kp is Fermi momentum 

in case of SNM, F(A:) = [ ^^+^^'^'+^^~^^'^' ] and J„ rep- 
resents the volume integrals of the isoscalar and the 
isovector parts of the M3Y interaction. The details of 
the present methodology may be obtained in Ref. [l6l| . 
However, for solving the Einstein's equations for stellar 
structure, we need to consider the total energy density 
including mass (also called, the mass-energy density) e 
which is related to the e and baryon number density p 
as e = p{e + m) where m (~ 938.919 MeV) is the aver- 
age of the neutron and proton masses in MeV unit. As 
the saturation energy per nucleon eg — —15.26 MeV and 
the saturation density pQ — 0.1533 fm^'^ are used in this 
work, the corresponding total energy density at satura- 
tion is £0=141.597 MeV/fm3=2.524 x lO^-* gm/cm^. Ob- 
viously, as these two parameters eo and po are extracted 
on the basis of information on finite nuclei, they put the 



2 



> 

(D 



200 



100 



DD-F 
KVOR 
DBHF 
This Work 



-50 




0.2 0.4 0.6 0.8 1 



> 

(D 



< 




P (fm" ) 



0.2 0.4 0.6 0.8 1 



3^ 



P (fm" ) 



FIG. 1. (Color online) (a) (Left) Density dependent behavior of symmetry energy (E^ 



p) determined by different 



approaches, (b) (Right) The EoS determined by different models (E/A vs. p) are shown for beta equilibrated charge neutral 
neutron star matter. 



two constraints to the high density behaviour of any nu- 
clear matter EoS. The pressure-density relationship for 
the present EoS is consistent with the experimental flow 
data confirming its high density behaviour. The val- 
ues of other important quantities, L, Ksi,m, Kasy and 
Kr, defined and calculated in Refs.[13, ^^^o agree 
extremely well with those extracted from experiments. 

The NSE given by Esymip) = 1) - e(PiO) has a 
value of 30.71 ± 0.26 MeV at the saturation density ob- 
tained from this calculation which satisfies one of the 
constraints on the high density EoS. At higher densi- 
ties the NSE (see Fig. 1(a)) using DDM3Y interaction 
peaks at p « 1.95po and becomes negative at p « 4.7po- 
In Fig.l the symmetry energy and energy per nucleon 
(E/A) are plotted as a functions of the baryon density. 
The equation of state for the /3-stable charge neutral neu- 
tron star matter (see Fig. 1(b)) is calculated numerically 
using Eq.(l) with /?-equilibrated proton fraction deter- 
mined from the NSE. These plots compare the symme- 
try energy functions and the EoSs determined by the 
present calculation and several relativistic models (e.g. 
DD-F, KVOR, DBHF) [?] for the neutron star matter. 
In Fig. 1(a), the NSE calculated by the phenomenological 
RMF models using density dependent masses and cou- 
pling constants (e.g. DD-F, KVOR) goes on increasing 
with density and never become negative. In the relativis- 
tic Dirac-Brueckner-Hartree-Fock (DBHF) approach, the 
NSE increases more rapidly with density indicating very 
large proton fraction at higher density. This shows an 
opposite trend to the NSE function determined from our 
EoS. The saturation density (po) used in DD-F, KVOR, 



DBHF and our EoS are 0.1469, 0.1600, 0.1810 and 0.1533 
fm~'^ respectively. So the DBHF uses considerably larger 
density whereas the po used by other two models are 
slightly different from the experimentally extracted value 
of 0.1533 fm~'^. The values of NSE at saturation den- 
sity calculated by DD-F, KVOR, DBHF and our EoS are 
31.6, 32.9, 34.4 and 30.71 MeV respectively. It is clear 
that DBHF slightly overestimates the value of NSE at po 
which is roughly around 30 MeV. 

A negative NSE at high densities implies that the pure 
neutron matter becomes the most stable state. Conse- 
quently, pure neutron matter exists near the core [2(| of 
the NSs. Although the present EoS is 'stiff' since the 
SNM incompressibility Koo = 274.7 ± 7.4 MeV, but the 
NSE is 'super-soft' because it increases initially with nu- 
cleonic density up to about two times the normal nuclear 
density and then decreases monotonically (hence 'soft') 
and becomes negative (hence 'super-soft') at higher den- 
sities [l^ [13] . This is consistent with the recent evidence 
for a soft NSE at suprasaturation densities [2l| and with 
the fact that the super-soft [22] nuclear symmetry energy 
is preferred by the FOPI/GSI experimental data on the 
TT^/ir^ ratio in relativistic heavy- ion reactions for the 
stability of NSs. 

The 13 equilibrium proton fraction xp [23*1 of a NS con- 
sisting of neutrons (n), protons (p) and electrons (e) is 
completely governed by the density dependent behavior 
of NSE. Contrary to the relativistic models like DD-F, 
KVOR, DBHF etc. this work does not support the fast 
cooling via direct nucleon URCA process as the maxi- 
mum x/j is 4.4% only. Recently it has been concluded 
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FIG. 2. (Color online) (a) (Left) The variation of NS masses (M) is shown with central density (pc). The results of this 
calculation are denoted as 'Static' and 'Rotating' for non-rotating and rotating NSs at the Keplerian limit, (b) (Right) Mass- 
Radius relationship is given for static and rotating stars at Keplerian frequency using this EoS. The other three plots (DD-F, 
KVOR, DBHF) present the same relationship for static star only. 



theoretically that an acceptable EoS of asymmetric nu- 
clear matter shall not allow the direct URCA process 
to occur in NSs with masses below 1.5 solar masses [1]. 
However, the possibility of fast cooling [l^, via di- 
rect hyperon URCA or any other processes enhancing 
neutrino emissivities like tt^ and K~ condensates may 
not be completely ruled out. Also a recent experimen- 
tal observation suggests [1^ high heat conductivity and 
enhanced core cooling process indicating the enhanced 
level of neutrino emission, which may be due to Cooper 
pairing. Further theoretical studies and sufficient obser- 
vational data are needed to shed some light on the cooling 
phenomenon of NS. 

Let us now explore the various properties of static and 
rotating NSs using the proposed EoS. To study the ro- 
tating stars the following metric is used: 

where the gravitational potentials 7, p, a and uj are func- 
tions of polar coordinates r and 9 only. The Einstein's 
field equations for the three potentials 7,_£, a have been 
solved using Green's function technique [27'-'29l| and the 
fourth potential uj has been determined from other po- 
tentials. All the physical quantities may then be deter- 
mined from these potentials [sO]- The matter inside the 
NS is approximated as a perfect fluid. Solution of the 
potentials and hence the calculation of physical quan- 
tities can be done numerically. The field equations for 



rotating stellar structure are solved by using our EoS fol- 
lowing the procedure adopted by Komatsu, Eriguchi and 
Hachisu |27l . [3g| . We choose the 'rns' code written by N. 
Stergioulas [3l| in calculating rotating as well as static 
NS properties. 

In Fig. 2(a), we have shown the mass of the stars as a 
function of central baryon density (pc). Our results are 
plotted for the static and Keplerian limit. This is obvious 
from the Fig. 2(a) that for the same mass comparatively 
less central density appears for the rotating stars due 
to centrifugal action. It may be noted that as angular 
frequency (ft) becomes greater, the structure of NS js^] 
gets changed not only because of centrifugal flattening, 
but also because it is taking place against the background 
of a radially dependent frame dragging frequency. For 
comparison we have also plotted the results of three other 
EoSs as mentioned earlier in the text. The maximum 
mass for the static case is about 1.92 Mq with radius 
~ 9.7 km and for the rotating case it is about 2.27 M0 
with radius ^ 13.1 km. So a mass higher than 1.92 Mq 
would rule out a static star as far as this EoS is concerned. 
The phenomenological RMF models DD-F and KVOR 
predict maximum mass around twice solar mass for non- 
rotating star. The relativistic DBHF model calculates 
the maximum mass ~ 2.33 Mq and therefore, the DBHF 
predicts massive NS even for static case. This is also clear 
from the Fig. 2(b) where the mass-radius relationships for 
all the above EoSs are shown. 

To summarise, we have presented a nuclear EoS at 
supersaturation densities which satisfies both the con- 
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straints from NS and heavy ion collision phenomenol- 
ogy. Our results show that with the stellar configuration, 
which contain a large fraction of /3 equilibrated NS matter 
with a thin crust is able to describe highly-massive com- 
pact stars, such as the one associated to the millisecond 
pulsars PSR B1516-h02B with a mass M=1.94l°-J^ M© 
(Icr) [si] and PSR J0751-H1807, with a mass M=2.1±0.2 
Mq to a 1(7 confidence level (and 2.ll;Q g Mq to a 2a 
confidence level) 0|. In the case of PSR J1748-2021B, 
a millisecond pulsar in the Globular Cluster NGC 6440, 
the measured mass is 1^=2.74;°^^ Mq (2ct) [s^. There 
are few other EoSs which can explain such a high mass 
for static case, however, they fail to explain the expected 
behavior of the NSE. We would like to mention at this 
stage that a star may not rotate as fast as Keplerian fre- 
quency due to r-mode instability [361] . There have been 
suggestions that r-mode may limit the time period to 
1.5 ms Ist']. However, pulsar rotating faster (e.g. PSR 
J17482446ad) than this limit has already been observed 
[38,j . Further observations and a better r-mode modelling 
may shed more light on this issue. 

Modern constraints from the mass and mass-radius- 
relation measurements require stiff EoS at high densi- 
ties, whereas flow data from heavy-ion collisions seem 



to disfavour too stiff behavior of the EoS. The data 
from massive NSs and pulsars may provide an impor- 
tant cross-check between high-density astrophysics and 
heavy-ion physics. The variation of pressure with density 
for the present EoS is consistent with the experimental 
flow data conflrming its high density behaviour. We flnd 
that the large values of gravitational masses (~2.0 Mq) 
for the NSs are possible with the present EoS with the 
SNM incompressibility K^o = 274.7 ± 7.4 MeV, which is 
rather 'stiff' enough at high densities to allow compact 
stars with large values of gravitational masses 2 Mq 
while the corresponding symmetry energy is 'super-soft' 
as preferred by FOPI/GSI experimental data. Thus the 
DDM3Y effective interaction which is found to provide 
unified description of elastic and inelastic scattering, var- 
ious radioactivites and nuclear matter properties, also 
provides excellent description of /3 equilibrated NS mat- 
ter to allow the recent observations of the massive com- 
pact stars. 
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